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Abstract 

We study a fiat brane soiution in an effective 5D action for cascading gravity and propose a 
mechanism to screen extrinsic curvature in the presence of a large tension on the brane. The 
screening mechanism leaves the bulk Riemann-flat, thus making it simpler to generalize large 
extra dimension dark energy models to higher codimensions. By studying an action with cubic 
interactions for the brane-bending scalar mode, we find that the perturbed action suffers from 
ghostlike instabilities for positive tension, whereas it can be made ghost-free for sufficiently small 
negative tension. 



I. INTRODUCTION 



The problem of cosmic acceleration and its possible explanation as a cosmological constant 
have led to a wide variety of models in theoretical physics. Higher-dimensional theories of 
dark energy, in which our Universe is viewed as a AD brane living in a higher-dimensional 
bulk, offer an interesting proposal towards understanding dark energy as a manifestation 
of the presence of extra dimensions of space-time. Much progress has been made in this 
field using the braneworld picture in which all standard model particles are confined to 
a 4D brane, while gravity is free to explore the bulk [lH3]. This makes it possible to 
have cosmologically large extra dimensions [2llll|5]. The Dvali-Gabadadze-Porrati (DGP) 
model [6], in particular, takes this idea to the extreme and considers our 4D Universe to 
be embedded in a 5D bulk of infinite extent. Despite being observationally disfavored [T]- 
[TU] , the normal branch of the DGP model is perturbatively ghost-free, in contrast to the 
self-accelerating branch pTHlG] . and thus represents a perturbatively consistent infrared 
modification of gravity in which the graviton has a soft mass. 

Infinitely large extra dimensions also offer a promising arena for realizing Rubakov and Sha- 
poshnikov's proposal [IT] for addressing the cosmological constant problem, namely that 
brane tension could curve the extra dimensions while leaving the 4D geometry flat. While 
tantalizing, this idea immediately fails if the extra dimensions are compactified, since 4D 
general relativity, and hence standard no-go arguments [IH], apply below the compactifica- 
tion scale. Moreover, obtaining a flat AD geometry with compact extra dimensions requires 
canceling the brane tension against other branes and/or bulk fluxes [19]. The situation 
is more promising if the extra dimensions have infinite volume. The weakening of gravity 
as it enters the higher-dimensional regime (combined with an intrinsic curvature term on 
the brane) at least suggests that vacuum energy, by virtue of being the longest-wavelength 
source, might only appear small because it is degravitated [20H22]- 

The generalization of large extra-dimension dark energy models to higher codimensions 
is important not only for the cosmological constant problem but also for their possible 
embedding into string theory [^ [^ . Previous attempts of such a generalization have been 
found to give rise to a divergent brane-to-brane propagator and ghost instabilities around 
fiat space [2^ 125]. Furthermore, for a static bulk, the geometry for codimension N > 2 has 
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a naked singularity at a finite distance from tlie brane, for arbitrarily small tension |20j . 

The cascading gravity framework [26VI3T] avoids these pathologies by embedding the 4D 
brane within a succession of higher- dimensional branes, each with their own intrinsic cur- 
vature term. The brane-to-brane propagator is regulated by the intrinsic curvature term of 
the higher-dimensional brane. Meanwhile, in the simplest codimension-2 case, consisting of 
a AD brane embedded in a 5D brane within a QD bulk, the ghost is cured by including a 
sufficiently large tension A on the (fiat) AD brane: 

A > '^mlMl , (1) 

where mg = M^/M^, and Mn denotes the Planck mass in D dimensions. This stability 
bound was first derived through the decoupling limit M5, Mg — )■ 00, keeping the strong- 
coupling scale Ag = (mgM|)^/'' fixed. In this limit, the 6D framework reduces to a local 
theory on the 5D brane, describing weak-field 5D gravity coupled to a self-interacting scalar 
field 71. The bound ([T]) was confirmed in |31j through a complete perturbation analysis in 
the full 6D set-up. 

The codimension-2 solution exhibits degravitation: the brane tension creates a deficit angle 
in the bulk, leaving the geometry fiat. Since the deficit angle must be less than 27r, the 
tension is bounded from above: 

A < 2ttM^ . (2) 

Since Mg is constrained phenomenologically to be less than ~meV, this upper bound is 
unfortunately comparable to the dark energy scale. Given its geometrical nature, however, 
this is likely an artifact of the codimension-2 case and is expected to be absent in higher 
codimensions. This motivated [SU] to study the codimension-3 case, consisting of a AD brane 
living on a 5D brane, itself embedded in a QD brane, together in a 7D bulk space-time. In 
the limit of small tension on the AD brane, such that the weak-field approximation is valid, 
[30] showed that the bulk geometry is non- singular everywhere (away from the brane) and 
asymptotically fiat, with the induced AD geometry also fiat. 

In a recent paper [32], we proposed a proxy theory for the full 6D cascading gravity model by 
covariantizing the 5D effective theory obtained through the decoupling limit. The resulting 
action is a 5D scalar-tensor theory, describing 5D gravity and the brane-bending scalar 
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mode (denoted by vr), coupled to a AD brane. The scalar field is of the conformal galileon 
type p5], with a cubic self-interaction term [TTl Since our brane is a codimension-1 
object in this case, the equations of motion are more tractable and allowed us in [32] to 
derive a rich cosmology on the brane. A similar strategy was used in earlier work |35j to 
construct an effective AD covariant theory, which was shown to faithfully reproduce much 
of the phenomenology of the full 5D DGP model. See [36VI39] for related work. 

The goal of this paper is to explore whether this effective framework also allows for flat 
brane solutions with tension and, if so, whether such degravitated solutions are stable. In 
particular, are the bounds Q and ^ reproduced in the effective theory? 

Remarkably, we find that our 5D theory allows for flat brane solutions for arbitrarily large 
tension, with the bulk geometry being non-singular. The cascading origin of the theory is 
essential to the viability of these solutions: if we let rriQ — )■ oo, corresponding to turning 
off the cubic scalar self-interaction, the bulk geometry develops a naked singularity a finite 
distance from the brane, as in |4Uj . 

Our mechanism for screening the brane cosmological constant relies crucially on n. In order 
for the theory to have a well-defined variational principle, the cubic self-interaction term 
requires appropriate interactions for tt on the brane, analogous to the Gibbons-Hawking- 
York term for gravity. In the presence of brane tension, these scalar boundary terms screen 
the tension, resulting in a fiat geometry. This is the interpretation of our mechanism in 
the Jordan frame, in which the scalar is non-minimally coupled to gravity. There is of 
course a similar intuitive explanation in the Einstein frame. There, based on the Israel 
junction conditions, one would expect that a large brane tension should imply large extrinsic 
curvature, and hence large {i.e. super-Planckian) bulk curvature near the brane. Instead, 
the scalar boundary terms effectively screen the tension, much like the screening of charges 
in a dielectric medium, resulting in a small source for bulk gravity. 

The screening mechanism we propose seems to resolve the problem with earlier self-tuning 
attempts. A perturbative analysis of this mechanism, however, shows that it is difficult 
to avoid ghosts in such a model for positive brane tension, while it is possible to obtain 
consistent ghost-free solutions for negative tension. We further find that the model is free 
of gradient instabilities, and scalar perturbations propagate sub-luminally along the extra 
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dimension. It is also worth mentioning that we only consider solutions in which the bulk is 
flat, hence we are working on a different branch of solutions than those studied in [32], and 
our results are in no way contradictory to 



We have organized our paper in the following way. After briefly reviewing cascading gravity 
in Sec. [Tl| we present the flat brane solution in Sec. Ill In Sec. IV we discuss perturbations to 



the screening solution around a flat background, and derive various conditions for stability, 
both in the bulk and on the brane. We summarize our results and discuss future research 
avenues in Sec. |Vl 

A comment on our notation: We use the mostly positive signature convention. Indices 
M, iV, ... run over 0,1,2,3,5 {i.e. the 4+ ID coordinates) and indices /i,z/, ... run over 
0,1,2,3 {i.e. the 3 + ID coordinates). We denote the fifth dimensional coordinate by 
y = x^. 



II. OVERVIEW OF CASCADING GRAVITY 



Consider a 6D cascading gravity model in which a 3-brane is embedded in a succession of 
higher- dimensional branes, each with its own Einstein-Hilbert action 

r4 r _]Vf3 



5'cascade = / d^XA/^— ^i?6 + / d^Xy/^^R^ 
^bulk J 4- brane 

+ [ d^a;^/^ (~^^^ + ^matter) , 

J3-brane V ^ / 



(3) 



where, as mentioned earlier, M^i denotes the Planck mass in D dimensions. The gravitational 
force law on the 3-brane "cascades" from 1/r^ to 1/r^ and from 1/r^ to 1/r^ as the Universe 
transitions from AD to 5D and ultimately to 6D at the crossover scales rrig ^ and rrig ^ 
respectively, where^ 

^5=^, -6 = ^. (4) 



^ Strictly speaking, the 40 — >■ 5D — )• 6D cascading behavior of the force law requires m^^ < rriQ^ , thereby 
allowing for an intermediate 5D regime. If m'^^ > ttIq^ , on the other hand, the scaling of the force law 
transitions directly from l/r^ to at the crossover scale ttiq^. 
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As mentioned in Sec. [T], this theory allows for degravitated solutions — a 3-brane with 
tension creates a deficit angle in the bulk while remaining fiat. Furthermore, the theory is 
perturbatively ghost-free provided the 3-brane tension is sufficiently large that ([T]) is satisfied. 

In the decoupling limit M5, Mg — )■ 00, with the strong-coupling scale 

Ae = imtMiy/' (5) 



held fixed, we can expand the action (|3j) around fiat space and integrate out the sixth 

dimension [32l The resulting action is local in 5D and describes weak-field gravity 

coupled to a scalar degree of freedom vr: 

Ml 
2 



C — ^' 5 / j5 

•J decouple ^ I u. X 



'bulk 

+ / d^X 

' brane 



1 27 



/If 2 1 



(6) 

where {Sh)MN = —^^huN /'^ + • • • is the linearized Einstein tensor. The scalar tt is the 
helicity-0 mode of the massive spin-2 graviton on the 4-brane and measures the extrinsic 
curvature of the 4-brane in the QD bulk space-time. An obvious advantage offered by the 
decoupling theory is that the 3-brane now represents a codimension-1 object, which greatly 
simplifies the analysis. On the other hand, its regime of validity is of course restrained to the 
weak-field limit and therefore of limited interest for obtaining cosmological or degravitated 
solutions. 

In [32] , we proposed a proxy theory for the full QD cascading gravity model by extending ([6]) 
to a fully covariant, non-linear theory of gravity in coupled to a 3-brane, 

27 

^6 



bulk 



lom 



+ 1 d'^aJi/^l ^i^4 + £matter)• 
>/brane V ^ / 



(7) 



This reduces to ([6]) in the weak-field limit provided that Vl{tx) 1 — 37r/2 for small vr. 
In |32], we chose Vl{ti) = e^'^'^l'^ and derived the induced cosmology on a moving 3-brane in 
static bulk space-time solutions. Interestingly, this choice corresponds in Einstein frame to 
the generalization of the cubic conformal galileon [33], whose structure is protected by 
symmetries. While the proposed covariantization of (|6]) is by no means unique, our hope is 
that ([7]) captures the salient features of the cascading gravity model, and furthermore 
that the resulting predictions are at least qualitatively robust to generalizations of ([T]). 



In this paper, we want to address whether ([T]) allows the 3-brane to have tension while 
remaining flat. To parallel the corresponding &D solutions, where the bulk acquires a deficit 
angle while remaining fiat, we will impose that the 5-D (Jordan-frame) metric is Minkowski 
space. For most of the analysis, we will leave Vtiji) as a general function, and derive con- 
straints on its form based on stability requirements. 

We work in the "half-picture", in which the brane is a boundary of the bulk space-time. 
In this case, the action ([t]) is not complete without the appropriate Gibbons- Hawking- York 
(GHY) terms on the brane [SI Il2], both for the metric and for vr |43j, to ensure a well- 
defined variational principle. These were derived in fiat space in ^] and around a general 
backgroud in [32], and the complete 5-D action is 

S =^ [ d'xV^,n (r, + K'- K,^K^^ + 2K^ - 2^ 

^ Jbulk \ \l \l 

6lm^ J bulk OZmg J brane V 

+ ^ / d^Xv^ -^i?4 + -Cmatter • (8) 
^ ^brano V ^ / 

Here q^^y = g^u—n^n^, is the AD induced metric, and K^j, = CnQfiuf^ is the extrinsic curvature 
of the brane, where n° is the unit normal to the brane, and Cn is the Lie derivative with 
respect to the normal. Note that we have added an extra factor of 1/2 in the brane action 
so that the Israel junction conditions obtained using ([s]) match with those obtained in the 
"full-picture" . The assumed Z2 symmetry across the brane guarantees that the bulk action 
in y > is equal to that in y < 0, while the bulk in (Isj) is defined only in y > 0. 
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lomg 



<9{A/(<9vr)^(9Ar)7r - -gMNdKid-nfd^'K - dMirdNnO^n 



Varying (|8|) with respect to the metric leads to the Einstein field equations 

1 
2 

- {gMND^ - VmVw) n , (9) 

where Gmn is the 5D Einstein tensor, and parentheses around indices denote symmetriza- 
tion: X(^MN) = {^mn + -^7Vm)/2. The matter stress-energy tensor on the brane is defined 
as 

r(4) ^ 2 6{^c 



Similarly, varying with respect to tt gives us the vr equation of motion, 

{D.nf - (VMdNTry - i?f ^Sa/vtS^tt = -^mln^R, , (11) 



--matter / 



with Q^T^ = dfl/dir. We further obtain the Israel junction conditions at the brane position 
by setting the boundary contributions to the variation of the action Q to zero. Variation 
with respect to the metric gives us the Israel junction condition 



(12) 



while varying with respect to vr yields the scalar field junction condition 

27 / \ 



(13) 



In the balance of this paper we seek flat brane solutions to the bulk equations ^ and (11), 
with boundary conditions set by (12) and (13). 



III. OBTAINING FLAT BRANE SOLUTIONS FOR ANY TENSION 



In this section we seek flat 3-brane solutions to the above equations of motion. To mimic 
the 6D situation where the brane remains flat but creates a deficit angle in a flat 6D bulk, 
we impose that the 5D (Jordan-frame) geometry is Minkowski space: 

dsbuik = VMNdx^^dx^ = -dr^ + dx^ + dy^ . (14) 

Similarly, the induced metric on the brane should also be flat. By Lorentz invariance, clearly 
we can assume the brane to be at fixed position, y = 0, with the extra dimension therefore 
extending from ?/ = to oo. By symmetry, we also have vr = 7r(?/). 

With these assumptions, the (5, 5) component of the field equations ^ and the vr equation 
of motion (11) are trivially satisfied, while the (/i, z/) components of ^ reduce to 

^" = J,r^ , (15) 



277r'2 
16m? 



n. 



where primes denote derivatives with respect to y. The junction conditions (12) and (13) 
can similarly be used to obtain the brane equations of motion. The vr junction condition 



(13) is trivial for a flat bulk and the (/i, z/) components of (12) reduce to, 

97r'3 A 
^ + ift^ = 2M| ■ 



where the subscript indicates that the function is evaluated at the brane position y = 0. 
We have further assumed that the matter energy-momentum tensor on the brane is a pure 
cosmological constant A, which we allow to be of any size, performing no fine-tuning like 
that usually required for the cosmological constant. In fact we would like A to be large (TeV 
scale), since we know from particle physics experiments that such energy densities exist on 
our 4D brane. Note that, although we neglect other matter for simplicity, its inclusion would 
not affect our overall conclusions. 



As a check, note that our junction condition ( 16 ) is consistent with the decoupling limit result 
ttq = A/3M| obtained in [26|,[3T]. Indeed, in this limit ^7 ,^0 ~ —3/2. Moreover, introducing 
the canonically normalized vr^ = M^^'^ir, we see that the ir'^ term drops out in the limit 



M5 — 7- 00, me — 7- keeping Ag = (mgM|)^/^ fixed. Hence our junction condition (16) 
reduces to the decoupling result in this limit. 



It is easily seen that the bulk equation (15) allows for a first integral of motion 



^ ttTt' H ^ = constant . (17) 

16mR 



Comparing against the junction condition (16) immediately fixes the integration constant 
in terms of A, and we obtain 

^ / 97r'3 A 



Notice that for suitable fi, (18) appears to admit a solution 7r(y) for arbitrarily large A. 



For example, suppose that A is large and positive, and we choose Q such that — )■ at 
large vr so that the cubic interaction term dominates everywhere, then this leads to a linear 
solution TT{y) increasing monotonically with y: 

Since vr is non-singular for any finite ?/, the solution is well-defined everywhere. Therefore a 
fiat brane solution is allowed for any tension. Of course, consistency of the effective theory 
requires that vr' <^ M5. Since vr is suppressed by the tiny scale mg, this is a weak requirement: 



TT 



Mr V 9M| ) \9mlMt 



« 1 , (20) 
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where in the last step we have used Q. Even with A ~ M4, this can be satisfied provided 
rriQ <^ m^. A hnearly growing 7c{y) is also desirable from the point of view of quantum 
corrections to the vr Lagrangian. It is well-known that such corrections are of the form 
(□tt)", that is, they always involve two derivatives per field, and hence vanish on a linear 
background. 

Note that the above remarks depend crucially on the cascading mechanism. If we let mg — >■ 
00, thereby effectively decoupling the sixth dimension and turning off the cubic vr terms 



in g, then ^ reduces to -Q' = A/2M|, with solution Q = -(A/2M|)?/ + c. For A > 0, 
as assumed above, the integration constant c must be positive since Q must always be 
positive (since it is the coefficient of i?5 in the action). Hence Q inevitably vanishes at some 
finite value of y in this case, indicating strong coupling. (In Einstein frame, this corresponds 
to a naked singularity.) The cascading mechanism, therefore, is crucial in obtaining a flat 
brane solution for positive tension. 

To gain further insight, we can translate to the Einstein frame: g^ff^f = MN- In this 
frame, the brane extrinsic curvature is non-zero and is determined by the Israel junction con- 
dition. Focusing on its trace for simplicity, and assuming fio = 1 without loss of generality, 
we have 



3 V IQml 2M| 



In the absence of the vr'^ term (corresponding to mg — )■ 00), the junction condition would 
imply K^/M^ ~ A/M^. In turn, requiring that the curvature remains sub-Planckian, <^ 
M5, would in turn impose a bound on the tension: A < [2U]. (Phenomenologically, M5 



must be less than ~ MeV, so this bound would be rather stringent.) Instead, using (16) 



and (20), we obtain 

M, \9mlMt ) ■ ^^^^ 

Again assuming mg <^ m^, this allows a Planck-scale tension, A ~ M4, while keeping 
<^ M5. In other words, the tt'^ contribution in (21) neutralizes the dangerous A term. 



leaving behind a much smaller curvature. This screening mechanism results in an effectively 
weak source for bulk gravity. This, however, also suggests that vr must be a source of negative 
energy to screen positive tension on the brane. This is not surprising since galileons are 
known to violate the usual energy conditions 
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Thus at the background level our proposed screening mechanism displays many desirable 
features. To be physically viable, the action ([T]) must be perturbatively stable around a flat 
bulk solution. We study this issue in detail in the next section. Unfortunately, we will find 



that the theory propagates ghosts around the large-tension solution (19). More generally, 
the absence of ghost instabilities, combined with the requirement that the bulk solution is 
well-defined everywhere, places stringent constraints on the form of Q and the allowed values 
of A that can be degravitated. In Sec. |V]we discuss possible ways to extend the framework 
to relax the stability constraints. 



IV. STABILITY 



In this section we study the stability of the degravitated solutions described above, by 
perturbing the complete Jordan frame action ([8| to quadratic order around the flat bulk 



metric (14). To do so, it is convenient to work in the Arnowitt-Deser-Misner (ADM) coor- 



dinates |l5] with y playing the role of a "time" variable, 

ds^5) = NW + Qt.Mx'' + N^dy){dx'' + N^dy) , (23) 

where A^ denotes as usual the lapse function and A"^ the shift vector. Focusing on scalar 
perturbations, we use the gauge freedom to make q^i, conformally fiat 

Moreover, we keep the brane at fixed position y = 0. (This of course does not completely 
fix the gauge in the bulk, but is sufficient for our purposes.) We perturb the lapse function, 
shift vector and scalar field respectively as 

N = 1 + 6N , (25) 
N, = d,f3 , (26) 
TT = 7i{y) + Tc{xf',y) . (27) 

Similarly, all functions of tt (such as ^{tt)) evaluated on the background will be denoted by a 



bar. (In particular, the background equations in Sec. Ill only hold for the barred quantities 
Tr{y) and ^(y).) 
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After some integration by parts, carefully keeping track of boundary terms, the complete 
action at quadratic order is given by 



pert 



Ml 



bulk 



6(^0' -6(^6N+ S^TT^ (dX) + + 8 



27M| 
32m2 



d^x 



6 -'bulk 

12 ca 



+ 27i"d'(3{n'6N - tt') 



+ 



Mi 



97M^ 



dS[27r'(a7r)' 



(28) 



'brane <JZ,//tg Jbrane 

Varying with respect to /3 and yields the first-order momentum and Hamiltonian con- 
straint equations, respectively. 



6N 
dX 



H TV 



^ 

TT TT ^ Vr Z 

TV Z vr 



where we have defined 



Z = ^ TT 

3 ' 



5m, 



2 • 



(29) 
(30) 

(31) 



Since 5N and j3 are Lagrange multipliers, either of the relations (29) and (30) can be 



substituted back into (|28|). The resulting quadratic action is 



S, 



Ml 



pert 



bulk 



+ 3M; 



d^x 



12fiC" + —{7^ + ^ 



brane 



Z2 V2f^ 



97r" 



5mJ 



(50^ 



-,-d\-—{W^7:) + —{dCy 



+ 



o7/\/f3 

d^a;[6(9C)'] - I dS[27r'(97r)2i 



brane 



32m? 



(32) 



6 J brane 



Note that the bulk action does not depend on vr, consistent with the fact that it is pure 
gauge from the bulk perspective. For consistency, its source at the brane position must 
vanish. That is, we must set the variation of the brane action with respect to vr to zero, 
thus obtaining 



TT 



1 - 



9 ^'^ 
Ami ^ 



(33) 
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Using this solution in (32) yields the complete C— action, 



S, 



Ml 



bulk 



- 3M| 



d^x 1 - 



Z2 
AmlZ 



^ 97r' 



97r'' 



Ml 



AmiZ^'"^' ■ 4 



d'x[6{dCy]. (34) 



'brane \ 'i/Zig^ / '±//igZv Jbrane 

where we have set = 1 on the brane, without loss of generality. As a check, we have 
repeated the bulk calculation in the Einstein frame, where the bulk geometry is warped, 
and obtained the same result. This calculation is presented in the Appendix. 



In order for bulk perturbations to be ghost-free, the coefficient of {dQ"^ must be negative: 

(35) 



Z^ Zn ^ Qvr" 

— + - < . 



2n 



n 



This inequality involves f2, vr' and vr". Using the background equations of motion (15) 



and (18), we can eliminate vr' and vr" in terms of and its derivatives, as well as the brane 



tension A. Hence (35) reduces to a second-order differential inequality for fi(vr), which 
constrains the allowed functions Q{7r) that can yield ghost-free solutions for a given value 



of A. More precisely, since (18) is a cubic equation for vr', we obtain up to three allowed 



differential inequalities for Q{Tf). The physically-allowed ^^(vr) should not only satisfy the 
ghost-free inequality, but must also be positive-definite and well-defined for all y > to 
avoid strong coupling. 

We have studied this problem numerically. Since it is non-trivial to solve the differential 
inequality directly, we have instead tried various forms for ^(vr) for different values of A, 



and checked whether these forms satisfied the ghost-free condition (35) for each of the roots 



of (18). For each root that satisfied (35), we then solved (18) for 7r(?/), and hence checked 
whether Cl{y) remained positive and well-defined everywhere. Some of the specific functional 
forms we have tried include = 1 ± 37r/2, e^^'^^'^ and 1 - 37r/2 + 97rV8. 

For positive tension, A > 0, we were unable to find any Q{tt) that could simultaneously 
satisfy the ghost-free condition and remain everywhere well-defined and positive. For large 



tension, A ^ Mg, any real root of (18) inevitably violates the ghost-free condition (35). For 
small tension, A ^ M|, it is possible to satisfy the ghost- free inequality, but the resulting 
fl{y) either vanishes or becomes cuspy a finite distance from the brane. This is illustrated 
in Fig. [l| for the case Q{7i) = 1 + 37r/2 and A = M|. 
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FIG. 1: In the left panel, we plot the quantity ^ + — |^ which appears in the ghost-free 



condition (35) for = 1 + 3tt/2 and A = M|. The three curves correspond to the three roots of 
the cubic equation (18) in vf'/me. The ghost-free condition requires ^ + — < 0, hence 



only the black (solid) curve is free of ghost instabilities. In the right panel, we plot ^{y) for the 
ghost-free case. Since O vanishes at finite y, corresponding to strong coupling, this solution is 
unphysical. We have found similar results for all positive values of A and functional forms of O 
that we have tried. 



For negative tension, A < 0, on the other hand, it is possible to find suitable ^(Tr) that 
satisfy the ghost-free condition and are well-defined for all y > 0. Figure [2] illustrates this 
for Q = 1 + 37r/2 and A = — Mg. However, this is only the case for sufficiently small values 
of the tension, |A| < Mg. For large values |A| ^ Mg, either the ghost-free condition cannot 
be satisfied or Q{y) is ill-behaved. The existence of non-singular, ghost-free degravitated 
solutions, albeit with negative tension, is certainly a welcome feature of our 5D covariant 
framework. That said, these solutions do not connect to the parent 6D cascading framework, 
where the deficit angle solution requires a positive tension source. 



Coming back to (34), there are other requirements that our degravitated solutions must 
satisfy. To avoid gradient instabilities in the extra dimension, the coefficient of ('"^ must be 
negative, which is automatically true since f2 > 0. Furthermore, from the ratio of the ('"^ 
and {d(y terms we can infer the sound speed of propagation in the bulk: 

= I zn.^ 9^ ' (3^) 



which is of course manifestly positive once (35) is satisfied. Using this we can determine 
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FIG. 2: Same as Fig. [ij except that Q, = 1 + 3it/2 and A = —Mq. From the right panel, we 
see that corresponding to the ghost-free branch is everywhere positive, hence this solution is 
physically viable. 



whether the propagation of perturbations is sub- or super-luminal. For the ghost-free ex- 
ample f2 = 1 + 37r/2 and A = — Mg shown in Fig. |2| is sub-luminal everywhere. 



Finally, the coefficient of {d()^ on the brane must be negative, in order to avoid ghost 
instabilities: 



-»-^0-SV'<O. (37) 



With 7715 > mg, for instance, this condition is satisfied for the negative-tension example of 
Fig. |2} As a check, we can compare this ghost-free condition with the stability bound ([T]) 
obtained both in the decoupling limit [26J and in the full 6D cascading framework [31j. In 
the decoupling limit with Q = 1 — 37r/2, where we expect agreement with the cascading 



results, (37) indeed reduces to A > 2mgM|/3. 



Note that the absence of ghosts on the brane can always be achieved by adding a suitably- 



large kinetic term for vr on the brane, thereby modifying (37) to a trivial condition. This 



intrinsic kinetic term would not affect the background solution nor the bulk perturbation 



analysis. In this sense, the bulk ghost-free condition (35) is a more robust constraint on the 
theory. 
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V. CONCLUSIONS 



Cascading gravity is an interesting approach to understanding dark energy as a manifestation 
of the presence of large extra dimensions. Unhke previous attempts, such as the DGP model, 
the propagators in cascading gravity are free of divergences, and the model has been found 
to be perturbatively ghost-free. Moreover, cascading gravity offers a promising arena for 
realizing degravitation: both in the codimension-2 |2nj and codimension-3 [20] cases, at 
least for small brane tension, the bulk geometry has been shown to be non-singular and 
asymptotically- flat, while the induced 4D geometry is fiat. 

In this paper, we have studied a recently-proposed effective 5D action of cascading gravity in 
an attempt to obtain fiat brane solutions. Our analysis has uncovered an intringuing screen- 
ing mechanism that can shield bulk gravity from a large tension on the brane, resulting in 
a small brane extrinsic curvature. The brane remains fiat for arbitrarily large tension, while 
the bulk is non-singular. Although this model offers an attractive mechanism to general- 
ize extra-dimension dark energy models to higher codimensions without any fine-tuning, the 
stability analysis imposes stringent constraints. The bulk solution is perturbatively unstable 
for positive brane tension, while it is possible to find stable solutions for sufficiently small 
negative brane tension. 

Our model agrees with earlier work in the weak-field limit, hence we do not contradict results 
that cascading gravity is indeed ghost-free. It does, however, raise the interesting question 
— is there a fundamental difference between a theory with large extra dimensions and an 
effective 4D scalar-tensor theory of gravity? A complete answer to this question demands a 
more detailed analysis, which we leave to future work. 

To improve stability, we are currently investigating the impact of including higher-order 
galileon terms for vr in the bulk, generalizing the results of [33] to 5D. Preliminary results 
show that these higher-order terms still allow for fiat brane solutions, while greatly alleviating 
the stability issues. In particular, ghost-free solutions are now possible with positive tension. 
However, demanding that gravity on the brane is approximately AD on sufficiently large 
scales appears to impose an upper bound on the brane tension. The results of this ongoing 
analysis will be presented in detail elsewhere. 
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Appendix: Alternative analysis of scalar perturbations 



In this appendix we present an alternative derivation of the bulk (^-action in (34), by per- 



forming the stability analysis in the Einstein frame: gf^j^j = Vl'^^^glfj^. We define a warp 
factor aE{y) = and a rescaled coordinate dyE = ^^^^dy. Removing the subscripts 

"i?" for simplicity, the bulk metric in Einstein frame is 



dsU^ = a\y){-dT' + dx') + dy' 



(Al) 



The Einstein frame bulk action is given by 

d^xy/^ 



, Ml 

•Jbulk — 



'bulk 



4 f^l 



2n 



~3n 



27Mi 



32m? 



6 -'bulk 



'71 



1 n f.^ 



(A2) 



Varying with respect to the metric yields the Einstein equations, M^Gmn = T^jj^, where 
the 71 stress-energy tensor, TJ^^ = —{2/y/—gr^)6S-„/6g^^^, is given by. 



rpTT 

-'-MN 



2M| n% 

3 fi2 



2dM7!-dNTl - gMN{d 



'7[ 



9M| fi,^ 
16ml iW^ 



gAmidn)^ - AdMTcdNiT{d7c) 



16m| 



d(^M{d7[fdN)-K - -gMNdK{d7rfd^7r - SA/vrSArTrDsTr 



(A3) 
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For the metric ( Al ) with -k = ^[y), the (5, 5) and (/i, u) components of the Einstein equations 
give us the following background evolution equations, 



6H^ 
3H' 



ip + p) 



(A4) 
(A5) 



where 



P 



TT ' 



12 



27 



5ms 



,1" /4 



27 



7r 



32m2' 



13 



n 



(A6) 
(A7) 



Here = a' /a is the Hubble parameter, with y playing the role of a "time" variable. 



To study scalar perturbations, we use ADM coordinates (23) and choose comoving gauge: 
<l^J.u = CL^iy)e^'^^^^'^^Vfj.L' ^ = 7r(y). In this gauge we cannot assume that the brane is at 
fixed position, but this is of no consequence here as we focus solely on bulk perturbations. 



The action (A2) can be rewritten using ADM variables as 



5i 



bulk 



Sq + S-ir 



(A8) 



with 



Mi 



Ml 



bulk 



bulk 



3 V / iV2 



27M| 



. TT 



12 



32m^ Jbuik 



3N 9 n 



where E^, = (g' - D^N, - D,N^)/2 = NK, 



(A9) 



Expanded to second order in the perturbations, 5N 
scalar field action reduces to 



M3 



bulk 



3H'^ + 3{4H^ + H') 



- 1 and = - 4:H, the 



-M\y)SN^ - M%y)SE''JN 



(AlO) 
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where 



27 



5m, 
27 



5m, 





-1/3 


6 






-l/3^'3 
' TC 


6 





"7r'^ + 7r'V + 12iJ7r 



3 Q 



/3 



(All) 
(A12) 



Varying the complete bulk action with respect to iV^ and N gives us the momentum and 
Hamiltonian constraint equations, 

' 2 



N 



/3 







(A13) 



^""^ - —H' + 3UH^ + H') + M^5N - M^SE'^o. = . (A14) 



For scalar perturbations, q^i, = a?{y)e'^^'^^^'y^r]^i, and A^"^ = d^(3, the first-order solutions 



to (A13) and (A14) are given by 

6C 



5N 

□4/3 



6 1 ^2 . -36i7' + AmM^ + AM^ - QM^ ^, 



(A15) 
(A16) 



As usual, we only need to solve the constraint equations at first-order in the perturbations 
to obtain the quadratic Lagrangian for (, since the second-order terms will multiply the 
unperturbed constraint equations, which vanish [16]. Also note that here d^P = q^^D^Dyfi 
whereas d'^C, = rf'^^d^dX- 



The quadratic action for ( is obtained by plugging back the solutions (A15) and (A16) into 



the original action (A8), (A9). We find that all of the □4/? terms add up to a total derivative, 



hence the final Einstein frame C~ action is 



Sr 



Ml 



bulk 



(A17) 



where 



B{y) 



Q{18H' - 2AHM^ - 2M^ + 3M^) 

{QH + M3)2 
Q{18H' + QHM^ + M^ + 3dyM^) 

{QH + M3)2 



(A18) 
(A19) 



and {dCf = V^'^d^CduC- 
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We can transform the action (A17) back to the Jordan frame by using the transformations 



between Einstein frame variables (now denoted with a subscript "i?") and Jordan frame 
variables: ue = fi^^^, dyE = i^^^^dy, and (e = C- The result is 



^Jordan 



Ml 



d^x 



bulk 



l2^]2 9it" 



Z2 \2n 



+ 



imi 



(A20) 



which matches with the bulk Jordan frame action in (34). 
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